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Abstract 


Extensive experimental research has indicated that active/inactive nanocomposites are promising electrode materials for rechargeable Li-ion 
batteries. Nanocomposite anode materials allow for capacities between 900 and 4000 mAh g~! whereas graphitic anodes, which are currently being 
used by industry, allow for a much lower capacity of 372 mAh g~!. By treating the active sites (which may be comprised of Si, Sn, Al, or Bi) as 
nanospheres embedded in an inert matrix, linear elastic fracture mechanics are employed in order to develop design criteria for these alternative 
battery systems, with respect to fracture that results from the large volume expansions that the active sites undergo upon Li-insertion. In particular, 
the present study: (i) predicts that smaller active site volume fractions are more stable; (ii) Griffith’s criterion is used to estimate the crack radius 
at which cracking will stop; (iii) based on the ultimate tensile strength of the inactive matrix the critical crack length at which the electrode will 
fracture is calculated; (iv) a theoretical estimation is made for the size of the active sites that will not allow cracks to develop and hence fracture 
of the electrode will be prevented. Based on the above analysis, Si active sites allow for a greater anode lifetime and therefore are preferred over 


Sn, furthermore the formulation can be applied to determine the most appropriate matrix materials. 


© 2007 Published by Elsevier B.V. 
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1. Introduction 


Lithium ion batteries are the main energy sources used in 
electronic devices such as cell phones and lap top computers, 
as well as in biomedical implantable devices. They are pre- 
ferred over other battery chemistries, such as nickel-cadmium 
and nickel-metal-hydride, not only because they are non-toxic, 
but also because their high energy density reduces their weight 
by half and their volume by 20-50%. Initially, both the cathode 
and anode were comprised of lithium metal but this resulted in 
safety issues, and therefore since the 1990s graphite is being 
solely used as the anode. Graphite, though, has a very low Li 
atomic density at full Li capacity in the carbon intercalation 
compound (LiCg¢), hence the volumetric Li capacity is rather 
low. This is an important aspect in choosing the electrode mate- 
rials because the total amount of Li-ions that can be taken up 
by the electrode material corresponds to the total time integra- 
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tion of the battery current. Therefore extensive research is being 
performed for the development of alternative anode materials 
that would give a higher energy density for the Li-ion battery 
systems. Some experimental research indicates that promising 
candidates for next generation Li-ion negative electrodes are 
engineered materials exhibiting a micro- or nano-scale compos- 
ite structure. This concept was originally suggested by Wang 
et al. [1] and later experimentally confirmed by Courtney and 
Dahn [2] and Thackeray et al. [3]. These nanocomposites are 
composed of materials capable of forming Li rich compounds 
(active materials) surrounded by a matrix material that is inert 
with respect to Li (inactive material). The matrix acts as a Li 
ion conductor while protecting the active sites from detrimental 
chemical reactions with the electrolyte. Both Si and Sn are good 
candidates for the active material in the composite structures 
as they have a very high Li intercalation density per host atom, 
forming the compounds Lig 4Si and Li4 4Sn. In terms of the total 
time integrated battery current per mass, these compounds have 
capacities of 900 and 4000 mAh g™!, respectively [4], whereas 
graphitic carbon gives a much lower Li capacity of 372 mAh g7! 
[5] Gt should be mentioned here that Si and Sn react with Li in a 
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similar way that C does). The drawback, however, is that during 
maximum Li insertion Si suffers a 310% volume increase [4]; 
in fact all metals that can act as active sites such as Sn, Bi, Al 
exhibit over a 100% volume increase after Li alloy formation. As 
a result significant internal stresses are produced and during the 
initial electrochemical cycles loss of mechanical integrity takes 
place at the active site/matrix interface; this leads to fracture of 
the electrode and a loss of electrical integrity [4]. 

From the aforementioned it can be seen that it is essential 
to model these alternative anode configurations from a mechan- 
ics point of view, before they can be used commercially. First 
steps towards this direction have been taken (i) by Aifantis 
and Hackney [6] who examined the elastic internal stresses the 
nanocomposite electrode materials undergo upon charging and 
discharging, and (ii) by Aifantis and Dempsey [7] who predicted 
the stability of the anodes upon maximum Li insertion. In partic- 
ular, in Ref. [7] various electrode configurations (disk-like active 
sites, fiber-like active sites, and spherical active sites) that are 
in the experimental stage were examined by employing linear 
elastic fracture mechanics. It was shown that spherical active 
sites allowed for the most stable behavior, and therefore only 
this configuration will be examined in the present work. Specif- 
ically, after determining the energy released during cracking as 
a function of the crack propagation distance, an expression is 
obtained for the stability index which allows the determination 
of the most stable active site volume fractions. Furthermore, 
Griffith’s criterion is used to estimate the crack radius at which 
crack growth will stop, while the critical crack length at which 
fracture of the electrode takes place will also be determined. 
Finally, by considering some additional physics, the size of the 
active sites that will result in no cracking will be predicted. 


2. Electrode configuration 


The configuration of the problem at hand is shown in Fig. 1. 
It should be noted here that the crumbling effect at the active 
site/matrix interface is significantly reduced for active sites that 
have nanometre dimensions, such as spherical nanocrystalline 
Si particles with a 10nm radius [5]. Additionally, this minia- 
turization is very effective from an electrochemical point of 
view since bulk Si reacts with Li at very large temperatures 
(400 K), whereas nanocrystalline Si can form Li alloys at room 
temperature [8]. 

It should also be noted that the capacity of bulk Si is signifi- 
cantly reduced after a few cycles [9], whereas at the nanorange 
it is closer to its theoretical value; i.e. 3000 mAh g7! have been 
achieved for 40 nm thin films, after 25 cycles [10]. Based on this 
experimental background the anode is taken to be comprised of 
Si nanospheres embedded in an inert, with respect to Li, matrix. 
The resulting analysis will provide criteria towards the selection 
of the most appropriate matrix material. 

As was mentioned in the introduction after the first few elec- 
trochemical cycles crumbling occurs at the Si/matrix interface. 
Since the matrix is more brittle from the Si it is assumed that it 
crumbles and forms a damage zone, as shown in Fig. 2. Since this 
region is severely damaged it is assumed to support only radial 
stresses (all other stress components vanish in that region), and 


Fig. 1. Electrode configuration: spherical active sites (shaded) of radius a are 
embedded in an inert matrix (blank) of radius b. This figure depicts two unit 
cells. 


Fig. 2. Configuration of unit cell used in analysis; a and b are the radii of the 
active site and matrix, A is the free expansion of the active site if it were not 
constrained by the matrix, p is the crack radius. 


can therefore be approximated by a number of radial cracks with 
length p-a. Furthermore, it should be noted that in Fig. 2, a and 
bare the radii of the Si and matrix, while A is the radial displace- 
ment to which the Si would expand to if it were not surrounded 
by the matrix. 


3. Energy released during crack growth 


The following formulation is a summary of the work in Ref. 
[7]. Upon electrochemical cycling an exchange of Li-ions takes 
place between the two electrodes. Upon charging of the battery 
Li-ions from the cathode are attracted to the negative electrode 
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and diffuse into the active sites, which as a result experience a 
300% expansion, while upon discharging Li-ions return to the 
positive electrode. Since this study is concerned with fracture of 
the anode the following analysis is valid once maximum charg- 
ing has been achieved. Upon maximum Li-insertion the stress 
that the Li-ions exert inside the active site is constant and there- 
fore the pressure that the active site exerts onto the matrix is 
also constant and is set equal to p. Similarly the pressure that is 
exerted onto the unit cell under examination by a neighbouring 
cell is constant and set equal to g. This implies the conditions: 


Orr(b) = —q a) 


where o,» is the radial stress. 
Since conditions of spherical symmetry are present the stress 
and displacement expressions in the active site can be written as 


1+ vs _ 2(1 + vs)r 
1 — 2v; ~ E 


where o,r, 099, Ogg are the stresses in the r-, 0- and -directions 
respectively, furthermore in Eq. (2) r is the radial distance from 
the center of the active site, u, is the displacement in the radial 
direction, vs denotes the active site Poisson ratio, E; the elastic 
modulus of the active site, and D, is an integration constant to 
be determined from the boundary conditions. It can be seen that 
there is no stress variation in the r-direction since the stress in 
the active sites is uniform. 

Considering Eqs. (1) with (2) allows the displacement inside 
the active site to be expressed as 


(1 — 2vs)rp 
E; ` 


Or(a) = =P, 


Orr = 060 = Oph = 2 Ds; uy Ds (2) 


r= (3) 
Now we can proceed to develop a displacement condition at 
the active site/matrix interface. If the active site was not sur- 
rounded by the matrix it would attain a radius of r=a+ A, upon 
maximum Li-insertion. For the present confined configuration, 
however, the surrounding matrix opposes this free expansion by 
pushing back the active site by a distance —ô. It can therefore be 
said that the initial radius of the active site is r; = a+ A, while the 
final radius of the active site, due to the presence of the matrix, is 
rg=a+ A — ô. Thus, the total displacement of the active particle 
surface, upon maximum Li-insertion, is Uur = rf — ri = —ô: 


(1 — 2us)(r + a)p 
Eg 


ur(a + A) = —ô = (4) 

Similarly the displacement at the active site/matrix inter- 
face can be formulated by considering the initial inner radius 
of the matrix annulus, wich is r=a and that which it attains 
after maximum Li-insertion which is r=a+ A — ô, according to 
the previous reasoning. Therefore at the active site/matrix inter- 
face (r=a)u, = A — 6. Combining this boundary condition with 
Eq. (4) gives 


(1 — 2vs)(A + a)p 
E; f 


ô= 


ur(a) = A (5) 

In the sequel the stress and displacement expressions inside 
the next region of Fig. 2, which is the damage zone will be 
developed. Since this region is severely damaged such that it 


can support only radial stresses, the stress equibrium relation 
gives 


do; 20/r 


k 
=0>o0,()=— 3, fora<r<p. (6) 
r 


dr r 
The integration constant, k, can be found by combining Eqs. 
(6) with (1), hence k= —pa’. Furthermore, from Eq. (6) the 

displacement expression is found as 
du, du, pa? pa? 


5 dr =e dr z Egr? t= Er 


Or(r) = E 


where Eg is the elastic modulus of the matrix. 

The constant of integration u“ is found by setting the dis- 
placement right in front of the crack tip u,() equal to u,(p), 
hence 

pa? : . par 
uCp) = Up) = — H u* = u* = uy (0) — (8) 
gP Egp 
where p is the distance from the centre of the active site to the 
crack tip, i.e. it is the crack radius. 

Insertion now of Eq. (8) into Eq. (7) gives the displacement 

expression inside the damage zone as 


pa? 


ur(r) = E 
g 


G z "| +uy(p), fora<r<p. (9) 
rp 

Finally, the displacement expression in the third region of 
Fig. 2, i.e. inside the uncracked matrix, can be found by treat- 
ing the matrix as a hollow sphere which experiences an internal 
pressure (by the active site) and an external pressure (by the 
neighboring unit cell). In fact, the solution for such a configura- 
tion has been solved by Westergard [11], and it reads 


r(1 + vg) | b3/(2r3) + (1 — 2ug)(1 + vg) 
Eg * b3/p3— 1 


u(r) = 


p? /(2r?) + (A — 2v,g)(1 + vg) 
1 — p3/b3 


\ foro<r<b. 
(10) 


The internal pressure p» is the pressure exerted at r=, and 
can therefore be found as p+ = pa?/p? (by direct substitution in 
Eq. (6)). 

It follows that the by letting r=p in Eq. (10) an analyti- 
cal expression can be obtained for u.(o), which can then be 
substituted in Eq. (9) for the development of a second bound- 
ary condition, in addition to Eq. (5), at the active site/matrix 
interface, as 


a> {1  21— vp) 
ua) = E : 
ae p 


3(1 — vgb? — Sbp?) 
2(b* — p*) 
(1) 
where $= qb?/(pa’). 
Finally, the displacement at the outer matrix boundary (r=b) 
is found by direct substitution in Eq. (10) as 


pa? | 3(1 — vgb? p — Sb) 
2bEg b3 — o 


u,(b) = +s+ wh -a3 
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Now we have all the parameters required in order to define 
the opening tensile stress responsible for crack growth and sta- 
bility. Since the damaged zone can support only radial stresses, 
and hence it consists of radial cracks, the stress responsible for 
cracking, is ogg. In fact, Dempsey et al. [12] have found that for 
an axial configuration such as the one at hand, ogg is given as 


pa? | 1 — 3S(p/b)? + 2(p/by? 


+) = 
b? | 2(p/by°( = (p/b)*) 


o(p (13) 


Based on Ref. [12] the energy released as the crack grows can 
be defined as 


2(1 — vpo% (pt) 


G(p) = ; 14 
(p) ni; (14) 
which in turn gives the stability index as 
b dG 
= ——. 15 
- G dp Ha) 


It can be seen that in order to define the crack opening stress it 
is necessary to know what p and S are. Since the battery system 
is under equilibrium, self equilibrated loading conditions imply 
that qb? = pa”, and hence $= 1. p can now be found by setting 
Eqs. (5) and (11) equal to each other, and solving for p, hence 


3(1 — vg) o? 
2(b? + bp + p”) 


a |p (1—vg) 
pE, |a 2 


p=A 


4 (a+ A)l — 2v; 


E, (16) 


The calculations performed in the sequel are concerned with 
this self-equilibrated loading case. 


4. Volume fraction considerations 


In manufacturing electrodes with the aforementioned struc- 
ture, it is of interest to examine the effect that the active site 
volume fraction has on crack initiation and propagation. In 
order to do so we must know the mechanical properties of the 
active sites. Upon Li-insertion in Si (and Sn) its crystal struc- 
ture changes from cubic to tetragonal, hence the elastic modulus 
depends on Li-content. Experimental data providing the modu- 
lus of active sites as a function of Li content are not available, 
however, so in the sequel the values of the modulus for bulk 
Si will be employed; this is a valid approximation if one con- 
siders that the purpose of the present work is to illustrate how 
to apply the above formulation to obtain design criteria for Li- 
batteries. Moreover, the ancillary parameters are not available 
at this time for many of the well-known lithium ion conductors 
that may act as matrix materials. Thus example matrix materials 
exhibiting a variety of mechanical behaviors will be exam- 
ined in order to develop specific design criteria for mechanical 
stability. 
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Fig. 3. Stability plots for systems with different active site volume fractions. 


4.1. Stability index 


In order to determine the size of the active sites that will 
result in a more stable battery system, we plot the stability index 
k with respect to the crack radius p, for b= 1 um and varying 
a. The various stability behaviors obtained are shown in Fig. 3. 
The active sites are taken to comprise of Si, while the matrix is 
of soda glass. 

The more negative the values at which the stability curves 
start out the more stable the system, since a greater energy dif- 
ference is required for crack growth to initiate. It can therefore 
be predicted from Fig. 3 that smaller active sites (i.e. smaller 
active site volume fractions f, where f= a?/b? result in more sta- 
ble anodes, since crack growth is initiated with greater difficulty 
for such systems. 


4.2. Griffith’s criterion 


Based on Griffith’s theory a crack will continue growing as 
long as the energy released (G) during its propagation is greater 
from the energy (Ge) required to create the new crack surface. 
By plotting therefore Eq. (14) together with the fracture energy 
(Gc) of the matrix material, which is a material constant, we can 
determine the crack radius at which the two energies intersect 
and hence estimate the distance at which crack growth will stop. 

In Fig. 4 we therefore plot G for various active site volume 
fractions together with Ge. The active sites are taken to be Si, 
whereas the matrix is Y203. It can be seen that the smaller 
the volume fraction of the active sites the smaller the distance 
at which crack growth stops. Hence, it is again predicted that 
smaller volume fractions of the active sites are more stable since 
they allow for less cracking. Furthermore, by plotting G and 
Ge, for the same volume fractions (i.e. a= 100nm, b=1 um), 
and active sites but different matrix materials, we can predict 
which matrix material allows for the smallest crack propaga- 
tion distance before cracking ceases; the results re shown in 
Table 1, where it can be seen that the most preferable material, 
based on these considerations, is Y2O3 (that is why it as used in 
Fig. 4). It is noted that in the construction of Fig. 4, the number 
of radial cracks, n, present had to be assumed, as Eq. (14) sug- 
gests. To obtain accurate values, respective experiments must be 
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Fig. 4. Griffith’s criterion for an electrode that comprises of Si active sites and 


Y203 matrix; various active site volume sites are considered. Arrows indicate 
the distance at which crack growth stops. 


Table 1 
Distance at which cracking stops for various matrix materials when a= 100 nm 
and b=1 pm 


Material Ge (m~?) Critical crack radius p at 
which cracking stops (nm) 

SrF2 0.36 725 

ThO2 2.5 450 

Y203 4.6 425 

KCl 0.14 840 


performed, since n, not only depends on the matrix material but 
it also changes according to p. Due to lack of such experimental 
data, however, and for illustration purposes, it is assumed that 
n=20. 

It should be noted that even if the actual Si elastic 
modulus value was employed in this section the qualitative 
patterns would remain the same, i.e. smaller active sites are 
more stable and hence crack growth initiation is more dif- 
ficult, and they also allow for smaller crack distances and 
Y203 is the more preferable matrix from those considered in 
Table 1. 


5. Critical crack length 
In addition to the above volume fraction considerations and 


use of Griffith’s criterion to determine when fracture will stop, 
direct use of the opening tensile stress Eq. (13), can provide 


Table 2 


information concerning the critical crack length at which frac- 
ture of the electrode will occur, as well as the size of the active 
sites that theoretically will result in no cracking. In particular, 
Eq. (13) gives us the tensile crack opening stress as a function of 
the crack radius (p), for any given material and geometric param- 
eters. If, however, we set ogg of Eq. (13) equal to the ultimate 
tensile strength of the matrix, and we define all the material and 
geometric parameters, the only unknown we are left with is p. 
Hence, we can solve numerically for p to obtain the critical crack 
radius and corresponding critical crack length (o — a) at which 
fracture of the electrode will take place for the corresponding 
system. 

In particular, it is of interest to keep the geometric parameters 
the same, while we vary the material parameters so that we can 
obtain which material selections allow for the greatest critical 
crack lengths. Such a sample comparison is done in Table 2; 
again Si is used as the active site, hence E£; = 112 GPa and v; = 
0.28, while various ceramics which are inert with respect to 
Li are used as the matrix. Two limiting cases are examined: (i) 
a= 10nmandb=1 wm, ie. the active site volume fraction is very 
small, Gi) a=750nm and b=1 um, i.e. the active site volume 
fraction is 0.42. It should be noted that the volume expansion of 
Si and Sn upon maximum Li insertion is approximately 300%; 
this corresponds to A =a! — a. 

Unfortunately, the same materials could not be used for both 
Tables 1 and 2 since for the materials in Table | the ultimate 
tensile strength could not be found in the literature, while for 
the materials in Table 2 the critical fracture energy (G,) could 
not be found. It can be seen from Table 2 that B4C is the most 
preferable matrix material since it allows for the greatest crit- 
ical crack radius and hence more electrochemical cycles can 
be performed before the electrode fractures. Furthermore, by 
comparing Sn and Si active sites it can be seen that the lat- 
ter are more preferable because they result in a longer anode 
lifetime since they allow for greater crack lengths prior to 
fracture. 

Furthermore, from Eq. (13) we can estimate the size of the 
active sites that will result in no cracking. When the crack radius 
(p) equals the radius of the active site (a) it implies that the 
crack length (po — a) is zero and hence no cracks are present. 
Therefore by defining the material parameters in Eq. (13), where 
again ogg = ultimate tensile strength, setting p =a and choosing 
a particular b we can estimate the corresponding radius a that 
will theoretically result in no cracking. Such sample results are 
shown in Table 3. 


Critical crack length for various ceramics (material parameters estimated from CRC Mater. Sci. Eng. Handbook [13]) 


Material Ultimate tensile Elastic modulus Poisson’s ratio Critical crack radius p (nm) for Critical crack radius p (nm) for 
strength (MPa) (GPa) active sites when a= 10nm active sites when a= 750 nm 
Si Sn Si Sn 
Al203 255 345 0.23 93 76 650.9 649.5 
B4C 155 450 0.21 123 99 651.6 650.9 
BeO 246 400 0.24 97 78 651 649.8 
WC 345 700 0.24 87 69 650.4 648.8 
ZrO 175 140 0.23 95 81 651.3 650.4 
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Table 3 
Active site size that results in no cracking when b= 1 um 


Material Radius a (in nm) of active site that 
results in no cracking when b= 1 wm 

Al0O3 757 

B4C 781 

BeO 746 

WC 745 

ZrO 758 


6. Conclusions 


The present study was able to develop design criteria 
(based on a theoretical mechanical analysis) for nanocomposite 
active/inactive anodes for rechargeable Li-ion batteries. Partic- 
ular conclusions are: (i) smaller active site volume fractions are 
preferable since they are more stable (i.e. crack growth initiation 
is more difficult), and they also allow for a smaller crack propa- 
gation distance to be attained before it is energetically favorable 
for crack growth to stop, based on Griffith’s criterion. This the- 
oretical prediction is consistent with the experimental evidence 
in Ref. [14], which shows that fibers comprised of nanospheri- 
cal SnO, crystals can be cycled without a loss in capacity. (ii) 
Si active sites result in a greater anode lifetime (compared to 
Sn) since their use allows for more electrochemical cycles to be 
performed prior to fracture. 

In general, it was illustrated that if the fracture energy (Gc) of 
the inactive material (i.e. matrix) is known, the distance at which 
cracking will stop can be obtained from Griffith’s criterion, while 
the critical crack length at which fracture will take place can 
be estimated, through Eq. (13) if the ultimate tensile strength 
of the matrix is known. Furthermore, for a particular materials 
system the size of the active sites that will theoretically result in 


no cracking can be estimated. Hopefully, the present work will 
motivate experimentalists to look further into the mechanical 
properties of materials being used as inactive matrices, and also 
measure the modulus of Si (and other active site materials) as a 
function of Li content. 
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